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ABSTRACT 

The objective of this project was to study the effect 
of small artificial cavities on natural convection from a 
horizontal surface. Tests were run with Feron 113. Data 
for heat flux as a function of bulk temperature difference 
were carefully obtained. These data yielded Nusselt number 
as a function of Grashof number or Rayleigh number. All of 
these were then compared with the data obtained by Hiep 
[Ref. 2]. Experimental results are presented for the heat 
transfer from horizontal circular disks, with and without 
arti f i ci al cavi ti es . 

The artificial cavities were found to affect the 
natural convection heat transfer from a horizontal surface: 
the heat transfer coefficients increased with the number of 
cavities. 
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NOMENCLATURE 



Symbol 


Descri pti on 


Units 


S 


Specific Heat 


BTU/lb -°F 
m 


D 


Test Surface Diameter 


Inches or ft. 


g 


Acceleration of Gravity 


Ft/sec 2 


h 


Convection Heat Transfer 
Coef f i ci ent 


BTU/Hr-Ft 2 -°F 


H 


Height of Cylindrical 
Tank 


Inches or ft. 


K ss 


Thermal Conductivity of 
304 Stainless Steel 


BTU/Hr-Ft-°F 


k f 


Thermal Conductivity of 
Freon 113 


BTU/Hr-Ft-° F 


P 


Power into Heater 


Watts 


R 


Radius of Cylindrical 
Tank 


Inches or ft. 


T b 


Bulk Temperature of 
.Fluid 


° F 


T s 


Test Section Surface 
T emperature 


° F 


6T 1 


Temperature Difference 
Between Test Surface 
and Fluid (T^-T^) 


° F 


T 2 


Temperature Measured in 
Stainless Steel Test 
Section about 0.72 
inch From Test Surface 


° F 


T 3 


Temperature Measured in 
Stainless Steel Test 
Section About 0.22 
Inch From Test Surface 


° F 
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Symbol 


Descri pti on 


Units 


6 T 2 


Temperature Difference in 
Test Section ( T £ - T ^ ) 


° F 


T f 


VV 

Film Temperature[ — ^ — ] 


°F 


6 X 


Distance in Stainless 
Steel Test Section 
Between Thermocouple 
No. 2 and Thermocouple 
No. 3 


Inches or ft. 


^subscript 


Uncertainly in variable 
denoted by subscript, 
for example: U T is 
1 b 

uncertainty in T^ 


Units of variable 
denoted by subscript 


a 


Fluid Coefficient of 


Ft^/sec 


3 


Fluid Coefficient of 
Thermal Expansion 


1 

° F 


y 


Fluid Dynamic Viscosity 


1 b m /hr-ft 


V 


Fluid Kinematic Viscosity 


Ft^/sec 


p 


Fluid Density 


lb/ft 3 

m 


Dimensionless 


Numbers 




Nu 


Nusselt Number Based on 
Test Section Surface 
Di ameter 


hD/k f 


Pr 


Fluid Prandtl Number 


y Cp v 
k f ’ a 


Gr 


Grashof Number Based on 
Test Section Surface 
Diameter 


g3 D 3 6T i 

2 

V 


Ra 


Rayleigh Number 


Gr. Pr 
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I. INTRODUCTION 



A. BACKGROUND 

It has been shown by Duncan [Ref. 1] that artificial 
cavities can have a significant effect on boiling heat 
transfer. If the cavities are of the correct size it is 
possible to induce nucleation at lower superheats than would 
otherwise be necessary. Recognizing the enhanced heat trans- 
fer which artificial cavities can induce in the boiling regime 
one is led to inquire as to the effect that these cavities 
might have on the natural convection heat transfer which 
precedes boiling. Hiep [Ref. 2] undertook an experimental 
investigation of the effect of artificial cavities on the 
natural convection heat transfer from a horizontal circular 
disc. Tests were run with water and Freon 113. For water 
O'Connor [Ref. 3] ran tests with a mirror finished test 
surface without the artificial cavities and plotted the data 
in terms of dimensionless parameters, Nusselt number versus 
Grashof number. Hiep started with seven artificial cavities 
for water, and gave the correlation result of Nusselt number 
in term of Rayleigh number as: 

Nu = 0.38 (Ra) 0,29 (1 ) 

For Freon 113 Hiep ran the tests with a blank and seven- 
hole test sections, again he gave the results in terms of 
Nusselt number versus Grashof number and Nusselt number against 
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Rayleigh number the correlation results for blank test section 
was 

Nu = 0.48 ( Ra ) 0 * 30 (2) 

The experimental apparatus that O'Connor [Ref. 3] and 
Hiep [Ref. 2] used was almost the same. Test surface had a 
diameter of 0.895 inches, installed in the center of a 
cylindrical tank 9.00 inches in diameter and 4.5 inches high 
and set at the same level as the tank floor (see Figure 3). 

A controllable resistance heater as a heat source was housed 
in a copper block attached to the bottom of test section. 

A calibrated resistor in series with the heater was used to 
determine accurately the input power to the heater. Eight 
thermocouples, wired in series were used to measure the 
average bulk temperature of the fluid and one removable 
thermocouple was used to provide a check on this temperature. 
To record the temperature distribution in the stainless steel 
cylinder, four stainless steel sheathed thermocouples were 
equally spaced along the axis of the test section. The 
temperatures at these locations were obtained using a Dana 
digital voltmeter. 

Since all the results and discussion of Hiep's work 
depend on the accuracy of determination of the test surface 
and fluid temperature the calibration of all thermocouples is 
extremely important. Hiep used the standard thermocouple 
tables instead of the average of his calibration value. The 
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uncertainty value in Hiep's thermocouple reading was 0.5°F, 
which results in a large uncertainty for heat transfer 
coefficient. The change in Nusselt number between the blank 
test section and the seven hole test section as determined by 
Kiep is within the experimental uncertainty assigned to the 
Nusselt number. 

The uncertainty in Hiep's experimental data, therefore, 
tends to make it difficult to draw a firm conclusion as to 
the effect which artificial cavities have on the natural 
convection heat transfer process. His experiments were 
limited to two test surfaces, blank and seven cavities. He 
also limited the range of surface temperature to those below 
the saturation temperature of Freon 113. In other words the 
experiment did not go into the boiling phase that Duncan 
[Ref. 1 ] obtai ned . 

B. THESIS OBJECTIVES 

The previous discussion has led to the objection of this 
experimental works, that is: 

1. To investigate the effects of the drilled cavities 
on natural convection heat transfer in Freon 113 
with greater number of test surfaces than Hiep 
obtained. 

2. To obtain results for the boiling phase of Freon 113 
to see if it agrees with the Duncan results. 

3. Try to keep the uncertainty value of heat transfer 
coefficient and other dimensionless parameters as 
low as possible. 



13 



II. EXPERIMENTAL APPARATUS 



A. DESIGN CONSIDERATION 

1. The experimental apparatus used in this study is 
the same as had been used by Hiep [Ref. 1]. A complete 
description can be found in his thesis and the thesis of 
O'Conner [Ref. 3]. A schematic diagram of the experimental 
set-up is shown in Figure 3. For convenience the ice-bath 
used by Hiep was replaced by an automatic electronic ice-point 
for the reference thermocoupl e. This automatic ice junction 
was set to 32°F as reference temperature. A Dana digital 
voltmeter that had been used to read and record the thermo- 
couples was replaced by a Hewlett Packard digital voltmeter. 

2. A Freon condenser was also added to recover the 
Freon that evaporates during the experiment. This condenser 
was connected to the top plate of the tank as shown in Figure 
2, and was cooled by tap water. 

B. TEST SECTION-HEATER ASSEMBLY 

Three test sections were used in this experimental work: 
one blank, one with seven holes which had been used by Hiep 
and one with thirteen holes (cylindrical drilled cavities). 

All three test surfaces had the same diameter, that is .896 
inch. Four 0.040 inch diameter holes were drilled along the 
length of each test piece to a depth of 0.375 inches, separated 
axially by about 0.5 inch and 90° apart from each other. 

These holes were for sheathed thermocouples installation in 
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Figure 



Photograph of Complete Apparatus 
Including Instrumentation. 
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Figure 2. Photograph of Assembled Experimental 
Unit. 
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Figure 3. Cross Sectional View of Tank, Test Section and 
Heater Assembly. 



Figure 4. Photograph of Blank Test Section 
After the Test. 
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Figure 5. Photograph of Thirteen-Hole Test 
Section After the Test. 
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Figure 6. Blank Test Section. 
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Figure 7. Seven-Hole Test Section. 
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Figure 8. Thirteen- Hole Test Section. 
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Figure 9. Heater Assembly. 
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each test section as shown in Figure 6-8. The thirteen-hole 
test section had thirteen cylindrical holes of the same size 
as the holes drilled on seven-hole test section, that is 
0.0375 inch in depth and 0.015 inch in diameter, drilled as 
shown in Figure 8. The heater section [Fig. 9] was the same 
as used by Hiep [Ref. 2]. 

C. INSTRUMENTATION 

1 . Temperature Measurement 

ISA type T copper-constantan thermocouples were used 
to measure the temperature distribution in stainless steel 
test section in order to get as accurate as possible in 
temperature distribution. These four thermocouples were 
sheathed in stainless steel (sheathed diameter = 0.035 inch) 
and the temperature of the surface was found by using a linear 
extrapolating method (see Appendix II, Fig. 18). 

The bulk fluid temperature was determined by using 
eight copper-constantan thermocouples wired in series (four 
pairs) manufactured using a Dynatech welder. These thermo- 
couples were wired in series into an insulated switching box 
so the reading gave an average temperature for four pairs of 
these thermocouples. A single removable thermocouple was 
used to provide a check on this temperature. The sheathed 
thermocouple leads were securely fastened to the bottom 
phenolic plate and then to a Jones bar arrangement. Because 
the sheathed thermocouples were difficult to manipulate and 
subject to twisting, and thus shorting, giving erroneous 
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reading, the more flexible wires were used for the bulk 
temperature thermocouples were used as extension wire into 
the switching box. A reference thermocouple was inserted 
into the automatic ice-point reference. A Hewlett Packard 
digital voltmeter accurate to 0.001 millivolt was used to 
read the output of the thermocouple in millivolt. 

2 . Power to Heater 

A Lambda regulated power supply was used to provide 
input voltage to the heater. A calibrated resistor was placed 
in series with the power supply and heater. This resistor 
was mounted on a large heat sink and was isothermal and 
therefore of constant resistance for all of the experimental 
runs. With this arrangement, a voltage reading across the 
resistor and across the heater were recorded. Knowing the 
voltage, and resistance, the current through the circuit could 
be determined by Ki rchoff current law: 



where 1^ = the current across the heater. 

= the current across the resistor. 
R = 2.013 ohm . 

= voltage across the resistor. 



(3) 



and the product of this current and voltage across the heater 
gave the input power, i.e., 

P ■ V H >H 
or 

P = V^(V^R) watts. (4) 
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This value could be accurately reset for different runs. 
The voltage across the heater and across the calibrated 
resistor were determined using the Hewlett-Packard digital 
voltmeter, with an accuracy of 0.001 volts. This power was 
not used to determined the heat flux but only to set the 
operating point from run to run. 



26 



III. EXPERIMENTAL PROCEDURE 



A. TEST SURFACE PREPARATION 

It was very important that the test surface be free from 
any machine marks and scratches that could influence the heat 
transfer other than the artificial holes. For this reason 
the surface was mirror finished. The machine marks were 
removed by using a Buehler Metallurgical Manual grinder; both 
0 emery polishing paper and 3/0 emery polishing paper were 
used. After that, the surface was washed with detergent and 
warm water to remove the previous abrasive and then washed 
with alcohol. In order to get a mirror finish the test 
surface was then wet polished on a metallurgical wheel covered 
with velvet. The wheel was impregnated with one micron diamond 
dust. The test surface was placed face down on the wheel and 
a low speed of rotation was used. During this work alcohol 
was used to wet the surface at all times. After completion 
of polishing, the test surface was washed with alcohol again 
and was dried with a dryer, and then wrapped with special 
paper to prevent contamination and touching. Before installa- 
tion the seven-hole and thirteen-hole test sections were put 
in an ultrasonic cleaner for ten minutes, were washed with 
distilled water, rinsed with acetone and dried. 

The seven-hole and thirteen-hole test surfaces were 
drilled before the mirror finish, by using a Sphinx spirec 
pivot drill, 0.015 inch diameter. The cavities depth were 



27 



controlled by using a feeler gauge on the lathe. A ratio of 
approximately 2.5 was desired for cavity depth to diameter, 
so the depth was about 0.040 inch; the location of seven 
holes and thirteen holes were as shown in Figures 7 and 8 
respecti vely . 

A height gauge which had a scale readable up to 0.0001 
inch was used to determine the location of the center of the 
thermocouple holes relative to the test surface. This gauge 
was also used to measure the diameter of the test surface. 

B. APPARATUS ASSEMBLY 

The main parts of the apparatus were the test section, 
heater assembly, and the tank assembly. The test section- 
heater assembly was securely bolted to the phenolic base 
plate, as described by Heip [Ref. 2] and once the insulation 
was in place and the can installed, the phenolic plate that 
served as the bottom of the tank was carefully fitted to the 
test section. The level of the test section was adjusted 
from the bottom base of the assembly by screws. The phenolic 
bottom plate of the tank was secured by three threaded rods 
120° apart, with the phenolic base plate. These three rods 
were also used to adjust the test surface in level with the 
phenolic bottom plate. For the tank assembly four "0" rings 
of 9.0 inch diameter were used for sealing to make sure that 
the cylindrical enclosure tank was leak proof. Two "0" rings 
for each side (top and bottom) of the tank were used. The 
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V 




cylindrical plexiglass container was then put in place and 
the top installed, six threaded rods, 60° apart were used to 
secure the top plate to the bottom plate of the tank. The 
glass laboratory condenser was installed on top of the test 
tank to recover any vaporized Freon. Two of the threaded 
rods were 18 inches long so they could be used to hold the 
condenser and insure that it would sit tightly on the 1 inch 
hole at the top plate of the tank. This 1 inch hole was used 
to pour the fluid into or out of the tank. To hold the 
condenser an adjustable holder was connected from each of 
the 18 inch rods to make the condenser stand still. The 
condenser was held tightly on the top but was also removable 
to allow filling of the tank. 

Six small holes were drilled through the plexiglass top 
of the tank to insert four pairs of thermocouples in the 
tank to measure the average bulk temperature of the fluid, 
for each pair one of the thermocouples was put near the 
bottom of the tank and the other one was placed near the top 
so that the average of fluid temperature all over the tank 
could be measured. Two holes, one in the middle of top plate 
and the other on the side of top plate, were used to insert 
one removable thermocouple to the middle of the tank to 
provide a check of the bulk temperature read by the other 
four pairs. 
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c. 



TEST PROCEDURE 



1 . T emperature 

A complete data run consisted of starting at the 
lowest heater setting, waiting until steady state was reached 
(usually about ten hours for each data point). The data was 
recorded every hour during the daytime, and at longer inter- 
vals during the night. After steady state was reached and 
the data for that power setting was recorded, the power to 
the heater was increased by approximately 5 watts. This 
procedure was repeated until boiling occurred from the surface. 
The power was then decreased in 5 watt increments and the data 
was recorded after steady state was reached. This procedure 
was repeated until the starting point was reached. The 
procedure of taking data during both a heating and cooling 
process was done to check on any hysteresis effects which 
might exist. 

The data recorded for each setting were the bulk 
temperature, the temperature of the individual sheathed 
thermocouples in the test s'ection, and the voltages across 
heater and calibrated resistor. When the bulk temperature 
and the temperature of the test section changed by less than 
half degree F per hour, steady state was assumed to be 
re ached . 
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2 . 



Determination of Nusselt, Grashof and Rayleigh 



Numbers 

For each run, at all input settings, the temperature 
distribution in the stainless steel test section are plotted 
and found to be linear. The surface temperature could be 
determined within 0.5°F from the plot of temperature versus 
thermocouple location by extrapolation, as shown in Appendix 
B. The properties of the fluid were determined at the film 
temperature which is the average of the bulk temperature 
(T^) and the test surface temperature (T<.) i.e., 



The thermal conductivity of 304 stainless steel 
was calculated using the data from Touloukian [Ref. 10] 
at the average of the temperature of the two nearest the rmo- 
couples to the test surface i.e.. 



length Sx. A simple energy balance at the surface using 
Fourier's Law of Cooling gives: 



T 



SS 



2 



This corresponds to the average temperature of the 




h6T 



h 



K SS * <5T 2 
<5 X • ( 6 T ) 



BTU/hr-ft 2 -°F 



(5) 
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The Nusselt number was obtained from the definition: 




( 6 ) 



Substituting values of h from (3) into (4) gives: 



Nu 




7Tx * ( 6 T -j ) • 



(7) 



The Grashof number was calculated using the diameter 
of the test surface as the characteristic length: 



Gr 



geo 3 d s -T b ) 

2 



( 8 ) 



The Rayleigh number is the product of the Grashof 
number and Prandlt number: 

Ra = Gr • Pr ( 9 ) 



where 




a_ 

v 




( 10 ) 
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IV. RESULTS AND DISCUSSION 



It was found that a minimum of ten hours was necessary 

for any input setting to reach steady state. Once steady 

state was reached, very good agreement between Nusselt number, 

Grashof number and Rayleigh number, based on test surface 

diameter, obtained during the heating and cooling phase of 

the run resulted. The uncertainty bands are shown on the plot 

of each curve. The investigation for the blank test section 

6 6 

covers a range of Rayleigh number from 30x10 to 350x10 , for 
seven-hole test section from 50x10^ to 250x10^ and 25x10^ to 
350x10 for thirteen-hole test section. A summary of the 
results for various data runs are shown in Tables I, II, and 
III. For each data run the Prandtl number, Nusselt number, 
Grashof number and Rayleigh number were calculated. These 
dimensionless parameters were determined using the test 
section surface diameter as the characteristic length. As 
described in the experimental procedure, data were taken at 
each setting as the apparatus was heated with increasing input 
and then again as the inputs were decreased, cooling the test 
surface. The cooling phase data are denoted in the table. 

Surface temperature of the test section was determined 
by linear extrapolation from all four sheathed thermocouple 
reading. This gave the surface temperature with an error of 
approximately .5°F. - 
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Summary of Results for Blank Test Section 
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Note: Represented the cooling phase. 



Summary of Results for Seven-Hole Test Section 



CD rd 
O OC 



to 

I 

o 

X 



CD S- 
O CO 



VO 

o 

X 

s- 

CO 



CD 3 

o 2 : 



LO 


CO 


CO 


00 






LO 


CD 


CXI 


r— 


00 


r— 


0 


CD 


(XI 


1 — 


CO 


r— 




0 


CXI 


00 


CO 


CO 


C\J 


O 






co 


co 


CO 


CO 


co 


CO 


CO 





CO 


.60 


.40 


p*~ 


LO 


CO 


CD 


.29 


LO 


« — 


co 




r— 


LO 


co 


CO 


CD 


1 — 


LO 


0 


VO 


O 




0 


VO 


O 


LO 




1 — 


r- 


CXI 


CXI 


CXI 


1 — 


1 — 





LO 


1 — 


1 — 


r^. 




vo 


1 — 


CD 


LO 


CD 


CO 


co 


CXI 


cxi 




co 


CO 


CD 


r^. 


r- 


CO 




LO 




CO 


r- 


r^. 


vo 
















vo 



CO 


0 


0 


0 


0 


0 


0 


0 


0 


CO 


cxi 


r— 


LO 


vo 


CXI 


0 


CO 




CXI 


LO 






vo 


CD 




r>v 


CXI 


vo 


CO 


1 — . 


vo 


CO 




, 


CO 


vo 




1 — 


CXI 


CXI 


CO 


CXI 


CXI 


I — 





LO 


0 

0 




CXI 


LO 


LO 


LO 


LO 




CXI 


CXI 


vo 


0 


vo 


1 — 




CO 




0 


r*“~ 


1 — 


CXI 


CO 


CXI 


1 — 


r- 


0 


CXI 


cxi 


CXI 


CXI 


CXI 


CXI 


CXI 


CXI 


CXI 



0 


LO 

vo 


vo 


LO 


vo 


r>* 


Cxi 


CD 


LO 


CO 


vo 




0 


vo 


0 






LO 


LO 


, 


vo 


CD 


t — 




CD 


vo 


0 


0 


CO 




LO 


CO 


vo 




CO 


1 — 


1 — 


1 — 


1 — - 


1 — 


CXI 


1 — 


r— * 


1 — • 


r— 



03 

CO 

rd 

JC 

CL 

CD 

C 



o 

o 



0 




CO 


0 


r— 




vo 


CD 


0 


0 




CD 




vo 


r— 






1 — 


CD 


<u 






00 


CO 


CD 


CO 


C 0 


CO 




SZ 



4 -> 

“O 

CU 

4 -> 

LO C 

COCOCOOOCOOLOO CU 



cn 


• 


• 
















to 


h- 


co 


0 


CO 


1 — 


0 


CXI 




CO 


0 


CU 




CO 


0 


1 — 


CXJ 


CO 


CXI 


r— 


0 


CD 








1 — 


r— 




1— 


r— 


1 — 


1 — 




CL 






















CU 






















u 






















* 














* 


* 


* 


* 




4-> 4-> 


CD 




1 — 






1 — 


vo 


CD 


LO 


• • 


0 4-> 


CXI 


0 


CD 


vo 






CD 


CD 


CXI 


CU 


CL rd 




CO 




0 


CD 


1 — 




CXJ 




+-> 


c 3 
















• 


• 


0 


►— ( S — - 


CXI 


LO 


r->H 


0 




0 




LO 


CXI 












r— 


f— - 


r— 











35 



Summary of Results for Thirteen-Hole Test Section 
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Note: * represented the cooling phase. 



A log-log plot of Nusselt number versus Grashof number 
and Rayleigh number based on test surface diameter is shown 
for each test surface in Figures 10 to 17. The curves show 
that the results fall into two regimes: 

1. Natural convection which is the purpose of this 
investigation. The Rayleigh number for this range falls in 

r r 

between 30x10 to 180x10 for these three test surface results. 
The curves show the straight line and almost have the same 
slope (see Figures 16 and 17). The correlation which results 
for the natural convection range of each test surface with 
Freon 113 as the fluid are: 

Nu = 0.48(Ra)^‘ 33 30x1 0^<Ra<l 80x1 0^ for blank test 

section. 

Nu = 0 . 53 ( Ra )° ’ 30 30x1 0 6 <Ra< 1 80x1 0 6 for seven-hole 

test section. 

Nu = 0.57(Ra)0' 3 ^ 30xl0^<Ra<180xl0^ for thirteen- 

hole test section. 

2. Boiling phase, the results of which must be 
considered with some caution. This is because the bulk 
temperature of the fluid was always below saturation and 
therefore the boiling washighly sub-cooled. In addition when 
nucleation occurred it began at the edge of the test surface 
in the space between the test section and the phenolic plate. 
Subsequently nucleation did take place on the test surface 

as can be seen by the cavitation damage visible in Figures 
4 and 5. The Rayleigh number for this range started from 
200x10 and up. The curve in this boiling phase is shown as 
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a, straight line with slope greater than the natural convection 
line (also see Figure 16 and 17). 

In the natural convection range the results show that the 
heat transfer coefficient of the fluid increases. This is 
consistant with Duncan's result but disagrees with Hiep's. 

Hi ep stated that as the number of drilled cavities increase 
the heat transfer coefficient decrease. The discrepancy with 
Hiep's result can be attributed to several factors: 

1. The uncertainty in thermocouple calibration which 
Hiep used was 0.5°F which led to a large uncertainty value in 
surrace temperature. Surface temperature was determined by 
linear extrapolation from two sheathed thermocouples near the 
surface. This has a large effect on the results which depend 
on the surface temperature. For the present experiment the 
value of the uncertainty in thermocouple calibration was 
reduced by carefully calibrating the thermocouples and using 
the average value of the calibration. This reduced the 
uncertainty in the thermocouple measurement to 0.2°F, which 
results in a smaller error in heat transfer coefficient. 

2. The number of test surfaces that Hiep used (blank 

and seven-hole) was insufficient to conclude that heat 

transfer coefficient decreases while the number of cavities 

increases. That is Hiep's results for the seven-hole test 

section fell within the uncertainty band of the blank test 

section. The present experimental work used a thirteen-hole 

test section in addition to the seven-hole test section. 

These results show a consistent trend of increasing Nusselt 
number with increasing number of cavities. 
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Log (Gr) 

Figure 10. Results for Blank Test Section Based on Grashcf Number. 
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Figure 11. Results for Blank Test Section 3aseci on Rayleigh Numbe 
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Log (fir) 

Figure 12. Results for Seven-Hole Test Section Based on Grashof flumbe 
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Figure 13. Results for Seven-Hole Test Section Based on Rayleigh Numbe 
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Log (Gr) 

Figure 14. Results for Thirteeen-Hole Test Section Based on Grashof Numbe 
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Fi ours 15. Results for Th i rteen -Ho 1 e Test Section Based on Rayleigh Numbs 
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Figure 16. Correlation of Results Based on Grashof .Number. 
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Figure 17. Correlation of Results Based on Rayleigh Number. 



V. CONCLUSION 



The results of this experimental work lead to the 
following conclusions: 

1. The artificial cylindrical drilled cavities affect 
the natural convection by increasing the heat transfer 
coefficient of fluid as shown in Figure 16 for the log-log 
plot of Nusselt number versus Grashof number and Figure 17 
for the log-log plot of Nusselt number versus Rayleigh number. 

2. The correlation results for this experimental work 
(the horizontal circular heated disk in cylindrical enclosure) 
are : 

o n 

Nu = 0.48(Ra)' , for blank test section. 

o n 

Nu = 0 . 5 3 ( R a ) ’ , for seven-hole test section. 

3 0 

Nu = 0.57(Ra)* , for th i rtee n - ho 1 e test section. 

All of these correlations were in the range of Rayleigh 
fi F\ 

number from 30x10 to 180x10 . For the value of Rayleigh 

f) F 

number from 200x10 to 350x10 the slope is increased as seen 
in Figures 10-17. 

3. The log-log plot of Nusselt number versus Grashof 

number and Nusselt number versus Rayleigh number shows that 

the curves were straight lines and had nearly the same slope 

in the range of Grashof number from 3x10^ to 20x10^ and 

6 6 

Rayleigh number from 30x10 to 180x10 . Between the range of 

Fi F) 

Rayleigh number from 180x10 to 200x10 the slope started to 
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change for each curve, and then the curves were straight lines 
again for Rayleigh number greater than 200xl0 6 as shown in 
Figures 16 and 17. 

4. Determining the surface temperature by the method 
of linear extrapolation as shown in Figure 18, gives the 
results within an error of 0.5°F. This is insufficient 
because the results show that if the surface temperature is 
changed by only .2-.3°F, the value of Nusselt number varies 
by about +_ 10% in the lower Dower input setting. This made 
the log of Nusselt number change by about +_! 2 % . 
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APPENDIX A 



THERMOCOUPLE CALIBRATION PROCEDURE 



The accurate determination of the test section and fluid 
temperatures was a major part of this experiment since all the 
results and discussion are based on these temperatures. For 
this reason, precise calibration of the sheathed and bulk 
temperature thermocouples was necessary. When the calibration 
was accomplished all the wiring and recording instruments that 
would be used in the experiment’ were not changed. 

A Rosemount Calibration System, with a constant tempera- 
ture oil bath, was used for the calibration. The ice junction 
thermocouple which was sealed in a glass tube filled with oil 
was placed in the electronic ice point reference junction. 

The four sheathed thermocouples and nine bulk temperature 
thermocouples were suspended several inches into the oil bath. 
A Platinum Resistance Thermometer in conjunction with a 
commutation bridge was used as a standard. The calibration 
was conducted over a range from 80°F to 220°F and a maximum 
error for the thermometer for this range of temperature was 
0. 004 mi 1 1 i vol t . 

The results of these thermocouples reading compared with 
standard thermocouple tables. The deviation from the table 
values ranged from +.004 to +.016 millivolt for sheathed 
thermocouple, and the deviation from table value ranged from 
+.006 to +.016 millivolt for bulk temperature thermocouples. 
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Since all the deviation is in the position range the decision 
was made to use the reading value and subtract the average 
difference from the table value. The values of uncertainty 
for sheathed thermocouples was 0.2°F and uncertainty of 0.2°F 
for the bulk temperature. 
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APPENDIX B 



SAMPLE CALCULATIONS 

For Thirteen-Hole Test Section: 

Input Power = 3.382 watts 
(Data recorded 6 October 1975) 

A. DETERMINATION OF THE SURFACE TEMPERATURE 

A graphical method was used to determine the surface 
temperature. For each input setting a plot of the temperature 
versus the thermocouple distance from the surface, as shown 
in Figure 18, was made using all four reading of sheathed 
thermocouples. The plots were linear and the surface tempera- 
ture could be obtained within an error of 0.5°F. 

B. CALCULATION OF THE PRANDLE NUMBER 
F rom the definition: 



Pr 



y. _ v 

a p 




( 10 ) 



where fluid properties were evaluated at film temperature. 

T c + T. 

[T f = -^-2 = 86 . 6° F] 

Fluid properties at T^ = 86.6°F [data from Table II from 
Thermophysical Properties of Refrigerants [Ref. 4]]. 
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T[°F] 

120 



Power Input 3.382 watts 
Bulk Temperature 79.4°F 
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90 I | l l I I I I I 

0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 

Distance from Test Surface, Inches. 

Figure 18. Determination of Surface Temperature. 
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V 



1.519 1 b^/ f t-hr 

k f = .04277 BTU/hr-ft- 0 F 

C n = .230 BTU/lb - ° F 

P m 



Using equation (10), we have 



Pr 



( 1 . 51 9 ) ( .230 ) 
.04277 



8.1686 



C. CALCULATION OF THE NUSSELT NUMBER 
Using the equations: 

_ K ss( 6l V , r , 

h " UxMfiT,) (5) 

and 

Nu = ^ (6) 

k f 

The value of K<.<^ was determined for an average tempera- 
ture of cylinder over the length 6x. The average of the two 
nearest temperatures to the test surface in the stainless 
steel cylinder were used. The value of K,~ s was obtained from 
a plot of temperature versus K ss based on values from 
Touloukian [Ref. 5]. Again fluid properties were evaluated 
at the film temperature using the data from Table II and 
Table 13 from Thermophysical Properties of Refrigerant 113 
[Ref. 4] and the following data were obtained for 3.382 watts. 
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105. 9°F 



T 3 = 100. 9°F 

K s$ = 8.926 BTU/hr-f t-°F 

D = .07467 ft. 

<ST 1 = 1 4 . 6 ° F 

6T 2 = 5 . 00° F 

6X = .04175 ft. 

P = 96.9 lb m /ft . 3 

The value of 3 were calculated by definition from the 
thermo-dynamic properties [Ref. 6] 



whe re 



1 dpi 
dT ' p 



1 / p 1 " p 0, 
p VV 



P 



P 



p 1 +p 0 
2 



T 1 - T q = Temperature difference of the given value 
p.j and p q . From Table 13 is 2°F. 



For this run , we have 



3 = 8. 668x1 0" 4 



Using equation (5) and (6), the Nusselt number is obtained: 



h 



8.926 ( 5. ) 
.04175(14.4) 



74.235 BTU/hr/ft 2 /°F 



74.235(0.07467) 

0.04277 



129.597 
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CALCULATION OF THE GRASHOF NUMBER 
Using equation: 

Gr = (gB/ v 2) D 3 6T ] 

The Grashof number is obtained 

Gr = 32.1 74 x 8.668x10 ( 0 . 07467 ) 3 ( 1 4 . 4 ) 

( 1.519 / 

^96.9x3600 

= 8.817X10 6 

CALCULATION OF THE RAYLEIGH NUMBER 
Using equation: 

Ra = Gr • Pr 



Thus 

Ra = (8.817xl0 6 ) (8.1686) 
= 7.202xl0 7 



APPENDIX C 



UNCERTAINTY ANALYSIS 

In this experimental work the method proposed by Kline 
and McClintock [Ref. 7] was used to obtain the uncertainties. 
Several assumptions and approximations had to be made in 
order to develop a reasonable estimate of the uncertainties 
for each variable. The following table is summary of these 
quantities. 



TABLE IV 

Uncertainty of Variables 



V a r i able 


Basis for Value 


Uncertainty (U) 


Cp 


assumption based on table 
values [Ref. 4J 


5x1 0~ 4 


BTU/lb m -°F 


D 


accuracy of measurement 


2x1 0~ 4 


ft . 


k f 


assumption based on table 
val ues [Ref. 4] 


1 xl 0” 3 


BTU/hr-f t-° F 


K SS 


assumption based on table 
values [Ref. 5] 


5 % 




T S 


based on scale of graphical 
method 


0 . 5° F 




T b 


thermocouple calibration 


0 . 2° F 




T 2 ,T 3 


thermocouple calibration 


0 . 2 ° F 




6 X 


accuracy of measurement 


4.2x10' 


4 ft. 


e 


assumption based on table 
values [Ref. 4J 


IxlO' 5 


1 / ° R 


p 


assumption based on table 
values [Ref. 4] 


lxlO" 2 


lb/ft. 3 

nr 


V 


assumption based on table 
values [Ref. 4] 


2x1 O -3 


lb /ft-hr 
m 
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The method that proposed by Kline and McClintock 
[Ref. 7] for determining the uncertainty values is the 
second-power equation. 

Uncertainties in the temperature differences 6T-| and 
6 T 2 were calculated from 

U 5 T, ■ v b = [( V 2 + <\ )2] * 

p 

= [ ( 0 . 5 ) 2 + (0.2) 2 ] 2 = . 538° F 

and 

U «T =U T -T = [dr + 0) T ) 2 ] % 

2 '2 '3 1 2 T 3 

2 2 * 

= [(0.2) + (0.2) ] = . 283° F 



Where U T = uncertainty in T 
1 S 



U T = uncertainty in T. 
'b b 



same as Uy U-j. = uncertainty in and respecti vely . 



Uncertainty in the convective heat transfer coefficient 
(h) was determined by using 



= [(• 



V SS 



‘SS 



) + ( 



6T ; 

TT7 



-) + (7^) 



6 X 



U 6 Tl 2 ' 



1 



( 11 ) 



All values of uncertainties are obtained from Table IV. 
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The Nusselt number is defined as 




( 12 ) 



and used of the second power equation for uncertainty in this 
number becomes 



2 % 



'N u 



-- C(-rr-) + (-#) + (r- 1 ) ] 



(13) 



It can be shown that the values of (-|~-) 
very small and can be neglected. Therefore, 



and 



are 





(14) 



To find the uncertainty of the Grashof number use this 
re 1 at i on : 

Gr = (gg/v 2 ) D 3 6T 1 



so 



Gr 

Gr 



U ; 

[</) 



+ <2^p) 



2 U 2 

+ ( 2 — ) 
V 



+ (3-D) 
D 



U 6T i 2 * 

] ) ] (15) 



6T 



1 



The values of U are found from Table IV. 
And by definition for Prandtl number is: 
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so 




U 2 U 2 U k 2 * 

[( T> + + <T7> ] 



M6) 



against all the values of U are found from Table IV. Find 

the uncertainty value for the Rayleigh number, 

where 



Ra = Gr *Pr 



Ra 

Ra 



- C(^r~) + (inf) 3 



Pr 



2 % 



Gr 



Pr 



( 17 ) 



A summary results of this analysis are given for each 
run in Table V for blank test section and in Table VI for 
seven-hole test section and for thirteen-hole in Table VII. 
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TABLE V 



Uncertainties in the Experiment (Blank) 



Power Input 
(Watts ) 


Nussel t 
Number 


Percent Uncertainty 
Grashof 
Number 


Rayl ei gh 
Number 


1.167 


20.09 


8.40 


8.70 


3.140 


9.12 


4.04 


4.67 


4.503 


7.33 


3.10 


3.90 


7.640 


5.90 


2.20 


3.30 


8.440 


5.80 


2.16 


3.23 


9.536 


5.69 


2.08 


3.19 


11.330 


5.55 


1 .99 


3.14 


24.500 


5.30 


1 . 89 


3.13 


29.655 


5.28 


1 . 83 


3.10 


34.735 


5.26 


1.80 


3.10 
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TABLE VI 



Uncertainties in the 
Power Input 



(Watts) 


Nussel t 
Number 


2.429 


11.30 


5.300 


6 . 80 


7.791 


5.98 


10.060 


5.66 


14.970 


5.42 



Experiment (Seven-Hole) 

Percent Uncertainty 

Grashof Rayleigh 



Number 


Number 


5.06 


5.64 


2.90 


3.74 


2.32 


3.33 


2.09 


3.19 


1 .95 


3.12 
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TABLE VII 



Uncertainties in the Experiment (Thirteen- Hole) 



Power Input 
(Watts) 


Nussel t 
Number 


1.167 


23.20 


3.382 


8.48 


5.758 


6.42 


9.176 


5.68 


16.300 


5.39 


26 . 280 


5.31 


34. 85 


5.29 



Percent Uncertainty 

Grashof Rayleigh 



Number 


Number 


10.26 


11.09 


4.00 


4.64 


2.75 


3.64 


2.15 


3.22 


1.98 


3.14 


1 .90 


3.13 


1 .91 


3.16 



62 



BIBLIOGRAPHY 



1. Duncan, D.S., "Natural Convection Heat Transfer From 
a Horizontal Disk in a Cylindrical Enclosure," M.S. 

Thesis, Naval Postgraduate School, Monterey, California, 
1971 . 

2. Le Vinh Hiep, "The Influence of Drilled Cavities On 
Natural Convection Heat Transfer From a Horizontal 
Surface," M.S. Thesis, Naval Postgraduate School, 

Monterey, California, June 1975. 

3. O'Conner, J.M., "Natural Convection Flow Visualization 
and Heat Transfer From a Horizontal Circular Disk," 

M.S. Thesis, Naval Postgraduate School, Monterey, 
California, June 1974. 

4. Thermophysical Properties of Refrigerants, 2nd ed., 

ASHRAE, 1973. 

5. Touloukian and DeHitt, Thermophysical Properties of 
Matter, V. 1 , 1970. 

6. Streeter, V.L., Fluid Mechanics, 5th ed., McGraw-Hill, 
1971. 

7. Kline, S.J. and McClintock, F.A., " Describing Uncertainty 
in Single Sample Experiments," Mechanical Engineering, 

V. 75, p. 3-8, January 1953. 

8. Marto, P.J., Moulson, J.A., Maynard, M.D., "Nucleate Pool 
Boiling of Nitrogen With Different Surface Conditions," 

J. of Heat Transfer, V. 90, p. 437-444, November 1968. 

9. Ostrach, S., " New Aspects of Natural Convection Heat 
Transfer," Trans. ASME , V. 75, p. 1 287-1 290 , 1953. 

10. Torrance, K.E., Orloff, L., and Rockett, J.A., 

"Experiments on Natural Convection in Enclosures with 
Localized Heating From Below," J. of Fluid Mechanics, 

V. 36, part 1, p. 21-23, 1969. 

11. Torrance, K.E., Orloff, L., and Rockett, J.A., "Numerical 
Study of Natural Convection in an Enclosure With 
Localized Heating From Below-Creeping Flow to the Onset 
of Laminar Instability," J. of Fluid Mechanics, V. 36, 
part 1, p. 33-54, 1969. 



63 



12. Holman, J.P., Heat Transfer, 3rd ed., McGraw-Hill, 

1 972. 

13. Knudsen, J.G. and Katz, D.L.," Fluid Dynamics and Heat 
Transfer," McGraw-Hill, 1958. 

14. Kays, VI. M., " Convection of Heat and Mass Transfer," 
McGraw-Hill , 1958. 



64 



INITIAL DISTRIBUTION LIST 



No. Copies 



1. Defense Documentation Center 2 

Cameron Station 

Alexandria, VA 22314 

2. Library, Code 0202 2 

Naval Postgraduate School 

Monterey, CA 93940 

3. Department Chairman, Code 59 2 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, CA 93940 

4. Assoc. Professor M. D. Kelleher 2 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, CA 93940 

5. First Lieutenent Apichart Penkitti, REA 2 



Military Research and Development Center 
Bangkok, Thailand 



65 






Thesis 

P3296 

c.l 



Thes i s 

P329& 

c.l 







* r \n o i 

k 'w JL 

Penki tti 

The influence of 
artificial cavities on 
natural convection 
heat transfer from a 
horizontal surface. 




Penki tti 

The influence of 
artificial cavities on 
natural convection 
heat transfer from a 
horizontal surface. 



thesP3296 

The influence of artificial cavities on 




3 2768 001 97954 5 
DUDLEY KNOX LIBRARY 



